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A scheme is discussed for measuring iVth-order coherences of two orthogonally polarized light 
fields in a single spatial mode at very limited experimental cost. To implement the scheme, the only 
measurements needed are the JVth-order intensity moments after the light beam has passed through 
two quarter-wave plates, one half-wave plate, and a polarizing beam splitter for specific settings 
of the wave plates. It is shown that this method can be applied for arbitrarily large N. A set of 
explicit values is given for the settings of the wave plates, constituting an optimal measurement of 
the TVth-order coherences for any N. For Fock states the method introduced here corresponds to 
a full state tomography. Applications of the scheme to systems other than polarization optics are 
discussed. 
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I. INTRODUCTION 

In recent years, interest in entangled states has grown 
significantly. One of the most successful experimental 
implementations has been achieved by entangling the po- 
larization degrees of freedom of twophotons in sponta- 
neous parametric down conversion Since then, a lot 
of effort has been devoted to producingentangled states 
with an ever higher photon number [2t6|- For the most 
part, entangled photon states are generated by postse- 
lection in such a way that every photon is found to have 
occupied an individual spatial mode Q. In that kind of 
setup, a quantum state tomography is usually conducted 
with the help of quarter- wave plates and half-wave plates 
which operate on every output port separately in such a 
way that each mode is analyzed along different orthogo- 
nal bases. The theory behind this method is intuitive and 
has been described exhaustively However, there exist 
experiments that generate polarization-entangled states 
of higher photon numbers in a single spatial mode 0, [loj . 
Here, it is possible to split the beam into as many spatial 
modes as there are photons to conduct a full state to- 
mography as in [H, but for higher photon numbers, this 
approach is experimentally very costly and inefficient. 

In this paper, we discuss a method that does not re- 
quire the beam to be separated in different spatial modes. 
It is based on a very general theorem, formulated by 
Mukunda and Jordan, which states that it is always 
possible to calculate the coherences of a photon field 
from photon correlation measurements in several differ- 
ent bases [ll|. This theorem has been used in proposals 
to measure all second-order coherences of a field, corre- 
sponding to the variances of the so-called quantum Stokes 
parameters [HI, 03 ■ Here we show that it is possible to 
measure all TVth-order coherences in a light beam con- 
sisting of two polarization modes with an arbitrary and 
unknown amount of photons in each mode. For this mea- 
surement we only require two quarter-wave plates, one 
half-wave plate, and a polarizing beam splitter, all acting 
on the same single spatial mode. The TVth-order inten- 



sity moment of that mode is measured after passa ge o f 
the photons through the mentioned optical elements [14| . 
If the incoming field is in a Fock state of N photons, 
this procedure corresponds to a full state tomography 
of that state 15[ . We note that the scheme is not lim- 



ited to polarization optics, but may also be applied to 
other two-mode systems where the necessary operations 
can be implemented, for ex ample p hotons in two different 
Laguerre-Gaussian modes [lg, LL2| • 



II. THE BASIS OF THE METHOD 

An often-used description for the polarization state of 
light are the Stokes parameters, which describe the state 
of a polarized light beam as a point on the Poincare 
sphere. For a classical coherent beam, the parameters 
are defined by: 
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where ot\ and a 2 represent the amplitude of the beam in 
two orthogonal linear polarizations. The quantum Stokes 
parameters are derived from the classical description by 
replacing the amplitudes and their complex conjugates 
with the annihilation and creation operators of the elec- 
tric field components. This leads to the following defini- 
tion for the quantum Stokes parameters: 
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As usual, for the electric field operators the bosonic com- 
mutation relation [aj,aj] = <% holds. The measurement 
of the quantum Stokes parameters involves the four av- 
erages (a\ai), (a\a 2 ), (a^ai), and (a\a 2 ). Schemes to 
measure these four quantities are standard textbook ma- 
terial [l8j , since the procedure is equivalent to determin- 
ing the polarization of a light beam and their values equal 
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those of their classical counterparts. However, the vari- 
ances of the quantum Stokes parameters, defined by the 
difference of two anticommutators 
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with i,j = {0,1,2,3}, are not equal to their classical 
counterparts. For example, for a coherent beam, in the 
classical picture all variances should vanish, whereas the 
analysis of the quantized fields shows that the variances 
actually obey an uncertainty relation and cannot all van- 
ish at the same time [l9| . They can be described in terms 
of all nine normally ordered second-order field correla- 
tions: (a\a\aiai), (a\a 1 aia2), {a\a\a\a\), {a\a\a\a2), 
(a\a\a20,2) , {a\a\a\a2) , {a\a\aiai) , and 
Out of those, only the three field corre- 
lations (a\a\aiai), {a\a\a\a2), and {a\a\a,2a2), which 
correspond to second-order intensity moment measure- 
ments, are directly accessible in experiments. 

Korolkova et al. first proposed a way to measure the 
diagonal variances Vu by using a half- wave plate and a 
quarter- wave plate [19| |. Later, in [l2j it was shown that 
it is possible to measure the variances of all quantum 
Stokes parameters (i.e., all second-order coherences of 
light) by conducting measurements of intensity-intensity 
correlations after the light has passed two quarter-wave 
plates and one half-wave plate for a set of nine specific 
positions of the three wave plates, a setup which con- 
stitutes a universal SU(2) gadget for polarized light and 
implements a general rotation in SU(2) space [20( - The 
knowledge of S and V already gives a good idea of the na- 
ture of the quantum state; however, the measurement of 
higher-order coherences may add even more information 
about that state. In particular, if the measured beam is 
in a photon-number state with N photons, that is, if the 
photonic state is of the form 



N 



\N- 



N 

E 



i, 



(4) 



the density matrix of that state has a size of (N + 1) x 



(N+ 1) and its (N + l) 2 elements correspond to all iVth- 
order coherences. Thus, for an iV-photon Fock state, the 
measurement of all iVth-order coherences is equivalent to 
a full state tomography. In the following, we show that 
with the setup depicted in Fig.[TJ which is slightly modi- 
fied with respect to the one discussed in [l2| , it is possible 
to determine all Nth-order coherences for arbitrary N by 
measuring only ./Vth-order intensity moments. 

The action of the SU(2) gadget on two orthogonally 
polarized modes a\ and a,2 can be parametrized by two 
angles 9 and </> 
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FIG. 1: (color online) A sketch of the setup with two quarter- 
wave plates (QPi and QP2), one half-wave plate (HP), and 
a polarizing beam splitter (PBS). At the detector (D), the 
iVth-order intensity moment is measured. 



where 9 and <f> are abstract parameters determined by the 
orientation of the three wave plates. The exact functional 
dependence is given in Scc. lIVI Using this parametriza- 
tion of the unitary transformation, we can express the 
most general case of measuring the .Vth-order correlation 
- defined by the correlation of the zth intensity moment 
in mode b\ with the (N — i)th intensity moment in mode 
62 - behind the SU(2) gadget as 



(bfb^-'bib^) = ]r J2 ( l ) ( ( N ~ *') ( N ~ 1 ) ( cos0 ) 2A 
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w J \y ) \ x j \ z 
x (sm9) w+x+v+z {-l) x+z e^ {x+y - w - z \a^ +x - w a\ N - M ^ (7) 



To solve for the (N + l) 2 independent real variables in 
the density matrix, we must perfom at least (N + l) 2 
measurements. Hereby, we must be sure that the val- 
ues of 9 and <f> chosen for these measurements lead to a 



system of independent linear equations; from Eq. Q, it 
is not obvious that this is possible for arbitrary N. In 
the following, a set of values for 9 and are be given for 
which we show that the measurement of ./Vth-order in- 
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tensity moments leads to a solvable system of equations. 
In the course of this proof, a natural recipe is developed 
that describes how the measurement results can be easily 
related to the coherences and populations of the initial 
state. 

III. A SET OF SOLUTIONS FOR EQ. © 

The results of this section show that it suffices to mea- 
sure photons of just one polarization, either b\ or &2, to 
determine all coherences. For this reason, it suffices to set 
up a measurement apparatus behind only one of the two 
output ports of the polarizing beam splitter (cf. Fig. [I]). 
In a random pick, we choose to measure the TVth-order in- 
tensity in mode b\ (i.e. i = N) and can therefore drop the 
summation over x and z in Eq. ([7]), which consequently 
simplifies to: 

(«)= f 0( N v )(^r- w - y 

w,y=0 ^ / \ » / 

(sin 9)™+y e ^y-^ (a\ N - w 4 w a^~ v a v 2 ). (8) 

Experimentally, (b\ b± ) corresponds to a measurement 
of the iVth-order intensity moment. Note that the num- 
ber of terms in Eq. © is (N + l) 2 and the expectation 
value of each coherence and population appears exactly 
once. Thus, the system of linear equations generated 
from this equation by (N + 1) 2 measurements of {b[ b^) 
for different <f> and 9 has exactly one solution if and only if 
we can choose the values of every pair (0, 9) such that all 
equations are independent. To arrive at such a choice, we 
first introduce new indices of summation, a and /3, such 
that we can rewrite Eq. © in a form where the phase 
e 1 ^^ factors out of one sum: 

(«)= e e (^V^Wr- 

0=-N a£Gi3 V 2 / V 2 / 

(smer(a{ N 2 4 2 af 2 a 2 2 ), (9) 

with 

a = y + w P = y — w. 



and 

G p = {2(N - k) - \P\} with k G {0, 1, . . . , N - \p\} 

(10) 

Equation (|9]) is the starting point for our analysis. Please 
note that for all k — 1 fcth roots of unity n (except unity 
itself), the equation X^=o( r K = holds. [25| This use- 
ful identity is exploited by choosing <p adequately to sim- 
plify Eq. ^ further and to introduce an inductive proof 
which shows that for suitable choices of 4> and 9, Eq. ([7]) 
can be solved. However, in order to arrive at a complete 
solution, we must distinguish in the following between 
measuring coherences of an odd or an even order N. 
A. N even 

If N is even, we choose for <j> the values <pk = with 
k e {1, 2, . . . ,N + 1}. Consequently, e ±l ^ k corresponds 
to all (N + l)th roots of unity. For every choice of (j>, we 
perform a measurement for N + 1 different values of 9, 
with 9j = and j £ {1, 2, . . . , N + 1}, thus carrying 

out (iV + 1) 2 measurements and obtaining (iV + 1) 2 differ- 
ent equations from Eq. (J9j> . By summing all equations of 
equal 9j , all terms from Eq. ([9]) with P ^ cancel because 
of the mentioned property of the roots of unity, and the 
sum over P contracts to P = 0. Thus, we are left with 
N + l equations (one for every value of j) containing only 
the N + l diagonal terms, each depending on a different 
power of cos 9j : 

{sm9 3 ) a (a\ N ~ § 4 f af~ f af), (11) 

where xg j = X^ fc x s- an< ^ x t k 1S * ne rcsu lt of the N 
photon measurement (b\ N b^) for setting 4> = <pk and 
9 = 9j. This set of equations can now be solved for the 
diagonal terms. 

We need not make more measurements to determine 
the other coherence terms. By first multiplying Eq. ([9]) 
by e l ^ k , and then adding all measurements for identical 
9j, only terms with p = —k and P' = N — k + 1 survive 
and we arrive at 
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Since N is even, all a are odd, while all a' are even or vice versa [cf. Eq. (fTU|) ] , leaving a total sum, in which 
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each coherence term again depends on a different power yields 
of cos 6j . [26] Furthermore, the total number of coherence 
terms appearing in Eq. (fl~2|) is given by \Gp \ + \Gpi | which 
is equal to N+l for every choice of k. Thus, the system of 
linear equations generated from Eq. (|12p by inserting all 



N+l values of 9j is solvable. Furthermore, (j)k determines 
the coherences that appear in the system. It is enough to 
generate a system of linear equations for every </>k with 
k G {1, 2, . . . , N/2} to solve for all coherences. Since 
the total number of measurements is equal to the total 
number of unknown variables, the presented set of values 
describes an optimal set of measurements. 



B. N odd 

If TV is odd, the previously described approach does not 
work, since a and a' in Eq. (|12p are both even or odd. 
Because of this, different coherence terms will depend 
on the same power of cos 9 and it is consequently only 
possible to solve for their sum. Therefore, we modify the 
choice of our values of <fi and 9 slightly: we choose A^ + 2 
settings for <f>, with <p k = k = {1, 2, . . . , N + 2}, so 



that e^ k describes all (N + 2)th roots of unity. For every 
4>k, we conduct measurements for N different values of 
9, with Oj = jy^-f and j G {1,2,..., N}. For this set 
of N(N + 2) measurements, we proceed as in the case 
for even N: the sum of all measurements for constant 9 
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which is identical to Eq. (fTTj) . However, we have only N 
equations to solve for A" + 1 terms, so we must conduct 
one more measurement [e.g. for (9, (f) = (0, 0)] to solve 
for all unknowns. At this point, the total number of mea- 
surements is again N(N + 2) + 1 = (N + l) 2 and, thus, 
also optimal. In the following, we need not make more 
measurements but can directly solve for the remaining 
unknown variables. In a first step, we multiply all equa- 
tions by e 1 ^ 1 before summation and arrive at 
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In contrast to the case for even N, we can arrive at a sys- 
tem of equations similar to Eq. (|13[) with only Af different 
terms, which we can solve immediately. Multiplying all 
equations with e l ^ k with 2 < k < (N+l)/2 gives all other 
necessary equations in a form equivalent to Eq. (|T2|) : 
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again with (3 = k, but j3' = N — k + 2. Since A^ is odd, 
all terms now depend on a different power of N, so that 
the system of linear equations again corresponds to a 
solvable (N + 1) x (N + 1) matrix, making it possible to 
determine all remaining coherences. For N — 1 (i.e., 
simply a polarization measurement), this leads to the 
choice of measuring the averages 



(o{oi), 

(a[ai) + {a\a 2 ) + e 1 ^ (a\a 2 ) + e~ l2 r {a\a\), 
(a\a{) + {a\a 2 ) + e ? ^ (a\a 2 ) + er 1 ^ (a 2 ai), 
(a\ai) + {a\a 2 ) + (a\a 2 ) + (ajai), 



IV. EXAMPLE AND MORE GENERAL 
APPLICATIONS 

In this section, we discuss the simplest nontrivial exam- 
ple, namely measuring all second-order coherences (case 
N = 2). For this task, we start by using the notation of 
Simon and Makunda [20( and write the unitary transfor- 
mation of Eq. ([5]) in terms of three Euler angles £, rj, £: 

U(£,r}>0 =exp (-i^Yj exp (*^) exp (~*"^r) ' 

" (17) 



which is different from what is discussed in standard text- 
books [lH, but equally optimal. 



where a 2 and 03 are Pauli matrices. Simon and Makunda 
show that the relation of the Euler angles to the actual 
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Euler angles (£, 77, £) 


(aQ Pl ,QQP 2 ,a H p) 


(iff, §ff) 
(iff, §ff) 

(Itt, fTT) 
(§ff,0) 
(|^,0) 
(!t,0) 


(1.775,0.676,4.197) 
(2.034, 1.318,5.176) 
(-3.833, 1.855, -0.692) 

(4.917,0.676, 1.775) 
(-1.107,1.318, -4.249) 
(5.591, 1.855,2.450) 
(0,0, -Itt) 
(0,0, -iff) 
(0,0, -fff) 


(1.673,2.011,0.169) 
(1.802,2.461,0.329) 
(2.010,2.938,0.464) 
(0.102,0.440, 1.740) 
(0.232,0.891, 1.900) 
(0.439, 1.367, 2.034) 

(3 71 "' i 71 "' if 71 ") 

(|tt, jtt, TgTr) 



TABLE I: All nine values of the Euler angles £, 77, and £ 
and the angles of the three wave plates in dependence of the 
parameters 8 and <j> f° r measuring second-order coherences. 



angles of the three birefringent plates is then given by [20| 
a QPl = (18a) 



«QP 2 
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with the Euler angles £, 77, and £ a function of the ab- 
stract angles and 



cos ^ = -\/a 2 + c 2 , 
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where a = Rc(e i *') sin0, b = Im(e 1 ' 1 ) sin0, and c = cos 0. 
In the case that a = c = 0, which occurs for 0, <fi = ±f , 
the corresponding Euler angles may be chosen as 77 = £ = 
and £ = 20, while in case that b = 0, which occurs for 
<p = 0,7r, the corresponding Euler angles may be chosen 
as 77 = £ = and C = -26* (C = 20) for = (0 = tt). 

With this translation of two abstract parameters into 
experimental quantities, it is now straightforward to cal- 
culate the settings for our wave plates for any arbitrary 
measurement. For example, if we wish to measure the 
nine variances of the Stokes parameters [Eq. <j3j) ] , the 
recipe from the previous section tells us to measure the 
second-order intensity after application of the nine uni- 
tary transformations that arise from all possible combi- 
nations of = iff, iff, |ff and (f> = |ff, |ff, 0. According 
to Eqs. (fl"8|) . every pair (0,0) corresponds to a certain 
triple of Euler angles (£, 77, £) and this in turn to a cer- 
tain triple of angles for the wave plates (QPi,QP2,HP), 
all of which are given in Table|H 

From this set of measurements of the second-order in- 
tensity moment, it is possible to calculate all variances 
of the Stokes parameters. In the case of a two-photon 



Fock state, this corresponds to all density matrix ele- 
ments (i.e., to a full state tomography). However, our 
method also serves for the determination of higher-order 
coherences of other (classical or nonclassical) states. For 
example, recently the covariance matrix of a Gaussian 
output state of an optical parametric oscillator has been 
measured plj |. If one would want to verify the Gaussian 
property of this state, the measurement of higher-order 
coherences like the ones discussed in the present work is 
required. 

The scheme presented here is also not limited to pho- 
tons of linear polarization: a\ and &2 may just as well 
correspond to any other pairwisc orthogonal photon po- 
larization modes. In fact, the general idea is applicable 
to any kind of bosonic multiqubit state where the equiva- 
lents to the needed devices exist: a universal SU(2) gad- 
get, a filter which transmits only one of the two qubit 
states, and a detector capable of performing a correla- 
tion measurement on the incident qubits. For example, 
one possible application could be the characterization 
of Laguerre-Gaussian beams with photons distributed 
among two different LG„ m modes [HI, H3]- I n this case, 
Agarwal discussed the SU(2) structure of their Poincare 
sphere [22l |: the equivalent of a polarizing beam splitter 
can be implemented with holograms [23j], and Ref. [24| 
points at the possibility of constructing an SU(2) gadget 
consisting of astigmatic lenses. 



V. CONCLUSION 

In conclusion, it was shown that a very simple exper- 
imental setup consisting of two quarter-wave plates, one 
half-wave plate, a polarizing beam splitter, and a mea- 
surement of higher-order intensity moments allows for an 
optimal measurement of arbitrary-order coherences be- 
tween two orthogonally polarized modes in a single light 
beam. Explicit formulas are given for the settings of the 
three involved wave plates. With these settings, the mea- 
surements allow the coherences to be obtained by a solv- 
able system of linear equations. The concept has been 
exemplified for the case N = 2, whereby, in the case of 
a Fock state, the capability of the method to perform 
a full state tomography has been demonstrated. The 
scheme could also be extended to include the measure- 
ment of phase-sensitive moments like (0102); however, in 
this case, one would need to add a local oscillator before 
the detector in Fig[TJ In a recent paper [2l|, the mea- 
surement of such phase-sensitive expectation values was 
reported for a Gaussian state. For the verification of the 
Gaussian property of such a state, the measurement of 
the A^th-order intensity moments like the ones presented 
in this paper is required. Finally, it was outlined that the 
method can be fruitfully applied to other systems as the 
general idea is not limited to linear polarization optics 
but is applicable to all bosonic systems where a univer- 
sal SU(2) gagdet and the analogs to a polarizing beam 
splitter and an intensity moment measurement can be 
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constructed. This property might make the present work 
interesting for a large range of similar topics. 



Acknowledgment 



U.S. thanks the Elite Network of Bavaria for financial 
support. 



P. G. Kwiat, K. Mattle, H. Weinfurter, A. Zeilinger, A. 
V. Sergienko, and Y. Shih, Phys. Rev. Lett. 75, 4337- 
4341 (1995). 

M. W. Mitchell, J. S. Lundeen, and A. M. Steinberg, 
Nature 429, 161 (2004). 

P. Walther, J.-W. Pan, M. Aspelmeyer, R. Ursin, S. Gas- 
paroni, and A. Zeilinger, Nature 429, 158 (2004). 

C. -Y. Lu, X.-Q. Zhou, O. Giihne, W.-B. Gao, J. Zhang, 
Z.-S. Yuan, A. Goebel, T. Yang, and J.-W. Pan, Nature 
Physics 3, 91 (2007). 

W. Wieczorek, R. Krischek, N. Kiesel, P. Michelberger, 

G. Toth, and H. Weinfurter, Phys. Rev. Lett. 103, 
020504 (2009). 

H. S. Eisenberg, G. Khoury, G. A. Durkin, C. Simon, and 

D. Bouwmeester, Phys. Rev. Lett. 93, 193901 (2004). 
N. Kiesel, Ch. Schmid, G. Toth, E. Solano, and H. We- 
infurter, Phys. Rev. Lett. 98, 063604 (2007). 

T.-C. Wei, J. B. Altepeter, D. Branning, P. M. Goldbart, 

D. F. V. James, E. Jeffrey, P. G. Kwiat, S. Mukhopad- 
hyay, and N. A. Peters, Phys. Rev. A 71, 032329 (2005). 
The two-photon case discussed in this paper corresponds 
to one photon in each of the two spatial modes. 

T. Tsegaye, J. Soderholm, M. Atatiire, A. Trifonov, G. 
Bjork, A. V. Sergienko, B. E. A. Saleh, and M. C. Teich, 
Phys. Rev. Lett. 85, 5013 (2000). 

L. K. Shalm, R. B. A. Adamson, and A. M. Steinberg, 
Nature 457, 67-70 (2009). 

N. Mukunda and T. F. Jordan, J. Math. Phys. 7, 849-853 
(1966). 

G. S. Agarwal and S. Chaturvedi, Jour. Mod. Opt. 50, 
711-716 (2003). 

A. Rivas and A. Luis, Phys. Rev. A 78, 043814 (2008). 

E. Shchukin and W. Vogel proposed a different technique 



[15 
[16 
[17 
[18 
[19 
[20 
[21 

[22 
[23 

[24 

[25 



to measure arbitrary-order field correlations of a light 

beam in a single polarization mode; E. Shchukin and W. 

Vogel, Phys. Rev. Lett. 96, 200403 (2006). 

R. B. A. Adamson, P. S. Turner, M. W. Mitchell, and A. 

M. Steinberg, Phys. Rev. A 78, 033832 (2008). 

M. Lassen, G. Leuchs, and U. L. Andersen, Phys. Rev. 

Lett. 102, 163602 (2009). 

M. T. L. Hsu, W. P. Bowen, and P. K. Lam, Phys. Rev. 
A 79, 043825 (2009). 

M. Born and E. Wolf, Principles of Optics (Cambridge 

University, Cambridge, 1999), 7th ed., Sec. 1.4. 

N. Korolkova, G. Leuchs, R. Loudon, T. C. Ralph, and 

C. Silberhorn, Phys. Rev. A 65, 052306 (2002). 

R. Simon and N. Mukunda, Phys. Lett. A 143, 165-159 

(1990). 

V. D'Auria, S. Fornaro, A. Porzio, S. Solimeno, S. Oli- 
vares, and M. G. A. Paris, Phys. Rev. Lett. 102, 020502 
(2009) 

G. S. Agarwal, J. Opt. Soc. Am. A 16, 2914-2916 (1999). 
A. Mair, A. Vaziri, G. Weihs, and A. Zeilinger, Nature 
412, 313 (2001). 

L. Allen, M. W. Beijersbergen, R. J. C. Spreeuw, and 
J. P. Woerdman, Phys. Rev. A 45, 8185 (1992). 
This can be seen instantaneously from the formula for 
the geometric series: 



fc-i 

E? 1 



[26] Remember that each coherence term appears only once 
in Eq. ((5} and thus also maximally once in Eq. (11211 , 



